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Abstract 

Motivated by recent inelastic neutron scattering experiments up to a high-energy region for 
La-based cuprates, we compute (q, u) maps of the imaginary part of the dynamical magnetic sus- 
ceptibility x(q> w ) i n the slave-boson mean-field approximation to the two-dimensional t-J model. 
While the strong spectral weight appears at incommensurate positions, namely at q ^ Q = (vr, 7r), 
for low energy, the incommensurate signals disperse with increasing uj and finally merge into a 
commensurate signal at a particular energy u = loq. These features are seen in both the <i-wave 
pairing state and the normal state. In particular, the incommensurate signals below loq in the 
normal state are due to the Fermi surface geometry, which we expect for La-based cuprates be- 
cause of a tendency to d-w&ve type Fermi surface deformations. Above cjq, strong signals appear 
to trace an upward dispersion especially for a low doping rate in the d-w&ve pairing state while 
typically broad spectral weight is obtained around q = Q in the normal state. Salient features 
of magnetic excitations in La-based cuprates are thus naturally captured in terms of particle-hole 
excitations. Global understanding of magnetic excitations in high-T c cuprates is discussed through 
a comparison with magnetic excitations in YBa2Cu30 y . 

PACS numbers: 74.25.Ha, 74.72.Dn, 74.20.Mn, 71.10.Fd 
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I. INTRODUCTION 



A major theoretical issue to understand the physics of high-temperature superconductors 
concerns their magnetic excitations. Most of inelastic neutron scattering measurements were 
performed for two high-T c cuprates, La 2 _. c Sr, c Cu04 (LSCO) and YBa 2 Cu 3 0j / (YBCO). For 
low temperature and low energy, incommensurate (IC) magnetic peaks are realized in both 



LSCO 



mperature ana low ene 
HQ and YBCO 0,0, 



], that is, the peaks of the imaginary part of the dynamical 
magnetic susceptibility x(q, appear at q = (tt, it ± 2^77) and (it ± 2ht], it); 77 parametrizes 
the degree of incommensurability. Above the superconducting transition temperature T c or 
possibly above the pseudogap temperature T*, the IC peaks merge into a commensurate 
H, _l while LSCO still retains IC signals up to a temperature much higher 



'0 3U 

Jab] 



peak in YBCO 
than T c . 

The size of the incommensurability 77 depends on the excitation energy to. For YBCO, 



r) decreases with increasing uj and vanishes at a specific energy u;q s . This commensurate 



signal was called a "resonance peak r 



11 



n 



12. 13] 

connected with the IC signals observed for lower energy, that is the peak disperses smoothly 



downwards to lower energy when q is shifted away from Q.|5j,|6 



14. 



151 ] Above the resonance 

energy, on the other hand, peaks of magnetic excitations trace an upward dispersion. 0, \i\ 
For LSCO, the energy dependence of 77 was 

a robustness of the IC signals was often contrasted with the behavior in YBCO. However, 
recent high energy neutron scattering data for LSCO with x = 0.10 and 0.16, 1^] and 
La 2 _ a ,.Ba :r Cu04 (LBCO) with x = O.I25Q revealed that the IC peaks disperse with uj. In 
particular, it was found in LBCO that the IC peaks merge into a commensurate peak around 

n 

55 meV,[19] above which upward dispersive features appear, similar to YBCO except that 
the data for YBCO were obtained below T c while those for LBCO above T c . 

Magnetic excitations in La-based cuprates are often interpreted in terms of the spin- 
charge stripe scenario. I20I] according to which (a tendency to) charge stripe ordering plays 
a central role. This scenario is based on the observation of two charge order satellite signals 
at both q = (0, ±4^77) and (±4^77, 0), whose wavevectors are just twice as large as those of 



20 



221 and 0.15. 



the magnetic IC signals, in Nd-doped LSCO with x = 0.10, |2l| 0.12. 
La9_T.Baj.CuO4 (LBCO) with x = 0. 125, (9J and Sr-doped LBCO with the hole density 



0.125. 



24L |25j It is assumed that each Cu0 2 plane has charge stripes characterized by a 
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single wavevector q = (0, ±4^77) or (±4^77, 0). Hence the observed two charge order signals 
are speculated to be superposition originating from different CuC>2 planes. Although such 
charge order signals were very weak and were observed neither for other doping rates nor in 
other cuprate superconductors such as YBCO and LSCO, the stripe scenario has attracted 
much interest . [26j 

There is another scenario based on conventional partic 
surface (FS), often refereed to as the fermiology scenario. 



27, 



role excitations around the Fermi 
MbaMbJ This scenario was 



explored in more detai 
tations in both YBCO 



and turned out to capture many important features of magnetic exci 



3y,flQflfly,y,HC, 42,0,0 and lsco.0,hq 



In particular, to be reconciled with a Fermi surface shape implied by angle-resolved pho- 
toemission spectroscopy (ARPES),0, 3 the fermiology scenario for LSCO was extended 



by invoking the idea of a <i-wave type Fermi surface deformation (c/FSD),|48( which is fre- 
quently called (i-wave Poo— in s taMity fl t fce FS e Xpmds alon g "e direct 
and shrinks along the k y direction, or vice versa. Comprehensive calculations of magnetic 
excitations for this d- wave-deformed Fsj^J showed that prominent features of magnetic 
excitations in LSCO are well-captured within the fermiology scenario. However, these cal- 
culations were confined to a low energy region of magnetic excitations. It is interesting to 
see whether the fermiology scenario can capture recently reported data for a high energy 
region in La-based cuprates also. Il9ll since such data were often discussed within spin-charge 
stripe— A^Q 

In this paper, we extend a previous calculation in Ref. 45], which was intended for mag- 
netic excitations in the low energy region for LSCO, to the high energy region. We compute 
the dynamical magnetic susceptibility x(q, w ) hi the slave-boson mean-field scheme of the 
two-dimensional (2D) t-J model. Combining the idea of the ciFSD with magnetic excitations, 
we show that salient features of magnetic excitations in La-based cuprates are well-captured 
up to high energy in terms of particle-hole excitations. We discuss how similarities and 
differences of magnetic excitations between LSCO and YBCO arise in the present theory. 

The paper is structured as follows. In Sec. II, we present the slave-boson mean-field 
scheme of the t-J model and show how the <iFSD is obtained in this scheme. The dynamical 
magnetic susceptibility x(q, uj) is formulated in the RPA with a renormalization factor. 



34, 



53J In Sec. Ill, numerical results of (q, to) maps of Inr\;(q, uj) are presented up to high energy 



for several choices of temperature T and hole density. In Sec. IV, we compare our results 
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with neutron scattering data in La-based cuprates and discuss similarities and differences 
of magnetic excitations between La-based cuprates and YBCO from the present theoretical 
viewpoint. Sec. V is a conclusion. 

II. MODEL AND FORMALISM 

We take the 2D t-J model on a square lattice 



defined in the Fock space with no doubly occupied sites. The operator cJ CT (c rcr ) creates 
(annihilates) an electron with spin a on site r, while S r is the spin operator. J(> 0) is a 
superexchange coupling between the nearest neighbor sites. We take into account hopping 
amplitudes up to the Zth nearest neighbors (I < 3), which are denoted by the conventional 
notation, t, t', and t", respectively. 

We introduce the slave particles, f ra and b r , as c rcr = blf ra , where f ra (b r ) is a fermion 
(boson) operator that carries spin a (charge e), and S r = \fl a cr apfrp with the Pauli ma- 
trices <r = (a x ,<j y ,a z ). The slave bosons and fermions are linked by the local constraint 
b r b r + fra fra — 1- This is an exact transformation known as the slave-boson formalism. 
We then decouple the interaction with the so-called resonating-valence-bond (RVB) mean 
fields:Q] XT=(Y,afrafva), and A T =(/ rT / r / i - f rl f r ^), with r = r' - r. These 

mean fields are assumed to be real constants independent of sites r. We approximate the 
bosons to condense at the bottom of the band, which leads to {b\.b r i) = 5, where S is the 
hole density. The resulting Hamiltonian reads 



with a global constraint J2 a {frafra) = 1 — 5; the k summation is over \k x ( y )\ < n. The 
RVB mean fields enter in £k and Ak, and are determined by minimizing the free energy. 
We obtain Xx = Xy = X an d A x . = — A y = A (d-wave paring). The explicit momentum 
dependence of the dispersion is given by 




(1) 





£k = —2 [t x cos k x + t y cos k y 
+ 2t'S cos k x cos k y + t"5 (cos 2k x + cos 2k y )] — fi , 



(3) 
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where t x = t y = t = t8 + | J% and // is the chemical potential, and the gap function is given 
by 



Ak = — - J A (cos k x — cos k y 



(4) 



The material dependence of high-7| cuprates is taken into account mainly by differ- 

0,0, 0, 



ent choices of band parameters. |2j 



cuprates. 



Following previous work for La-based 
29L l45j we choose tj J = 4, £' jt = —1/6, and t"/t = 0. In the course of the present 



work, we noticed that incommensurate structures of Im%(q, u) in the normal state become 
clearer for a larger t/J, showing better agreement with neutron scattering data. Note that 
t he Pres e, lt choice of t/ J is — with ao ab m U,o calcolatiooQ , which showed that 
realistic values for t/J lie in the range 2 — 5. 

The obtained FS is electronlike (holelike) for 5 > 0.10 (5 < 0.10) [Fig. Ufa)]. This FS 
seems to be inconsistent with ARPES data where a holelike FS was implied in a much 

n 

wider hole-doping region (6 < 0.20); |46j the reported FS looked nearly the same as a FS in 
YBCO.j3| If this were the case, a fermiology scenario predicts that magnetic excitations 
in LSCO should be essentially the same as those in YBCO, which contradicts with neutron 
scattering data especially for the normal state. This problem was considered in the 2D 
t-J model. It was found that the system has a tendency to a ti-wave type Fermi surface 
symmetry breaking. [48] The same tendency was also found in the Hubbard model. This 
CcFSD tendency is generated by forward scattering interactions of electrons close to the FS 
around (it, 0) and (0, it), and may be a generic feature as seen in various works: the slave- 

nn n 

boson mean-field theory,|48, 611 the exact diaeonalization analysis, 1 621 and the variational 

n nnn 

Monte Carlo study[63] of the t-J model, various renormalization group studies|6JJ, IfMj, |6Ja] 



and the renormalized perturbation the ory j6fil| of the Hubbard model, and the mean-field 
theory of the extended Hubbard model. |67l| 

The ciFSD competes with d-wave superconductivity. Re- 



model, 



48. 



61 



60 



61 



63, 



64, 



65 



66( In the t-J 



631 ] the d-wave pairing instability is dominant and the spontaneous (iFSD does 



not happen. However, the system still retains appreciable correlations of the dFSD,[48, 



n 



61 



l63| which make the system very sensitive to a small external anisotropy, leading to a 
strongly deformed FS. A natural origin of this anisotropy lies in crystal structures. La-based 
cuprates have the low-temperature tetragonal lattice (LTT) structure. It yields a 

small xy anisotropy, the direction of which alternates along the z-axis. Through a coupling 
to the LTT, therefore, we can expect an alternate stacking of a strongly deformed FS as 



shown in Fig. El Because of a weak interlayer coupling, the resulting Fermi surfaces have two 
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sheets, an inner electronlike FS and an outer holelike FS. The FS reported by Ino et al. 
may correspond to the outer holelike FS; the inner electronlike FS has not been detected 
successfully by ARPES. 

Although the deformed FS can be determined in a self-consistent calculation, we 
here deform a FS by hand by introducing a parameter a 

t x ( y ) = cd (5) 

with keeping t y r x ) = t. For simplicity we do not consider an interlayer coupling and compute 
magnetic excitations for the superimposed FSs shown in Fig. [Hb). We first determine 
RVB mean fields self-consistently for a = 1, and then tune a to get a FS compatible with 
ARPES; 46] we choose a = 0.85 for 5 = 0.15 [Fig. ^b)]. Since we here tune a by hand, the 
resulting hole density slightly deviates from 5 = 0.15 to be 0.156. 

We investigate magnetic excitations for both FSs shown in Figs. E^a) and (b). The 
irreducible dynamical magnetic susceptibility Xo(q ; ^) is given by 

E\r . E\r*r, \ 1 



Ynfq, uj) = — rr Ckk^ ( tanh — ^ — tanh k J" q ) , 
+ \ c k k+ q ( tanh ^ + tanh ^r^j 



x — - — ■ — + 



£ k + £ k+q + ^ + ir E k + E k+q - lo - ir 

where £ k = vf+A^, T is a positive infinitesimal, and 



(6) 



In a renormalized random phase approximation (RPA) the dynamical magnetic sus- 

ceptibility x(q, is given by 

/ \ xo(q, 
x(q, ^) = i , t/ \ — / — \ , s 
i + ^(q)xo(q, w) 

where 

J(q) = 2rJ(cosq x + cosq y ) (9) 

with a renormalization factor r. In the plain RPA one has r = 1, which leads to a divergence 
of x(q, 0) around q ~ (vr,7r) in a wide hole-doping region (5 < 0.17) in the ci-wave pairing 



state, signaling an instability to the antiferromagnetic (AF) state. However, fluctuations not 
included in the RPA obviously suppress the AF instability as shown in several numerical 



studies. 



• nl 22, Q 



This aspect may be roughly taken into account in a phenomenological 



way by setting r < 1. 



l53| Here we choose r = 0.35, which confines the AF instability to 



5 < 0.02, consistent with the phase diagram of LSCO.|74j| 
III. RESULTS 

We compute xi^ a; ) numerically for the FS shown in Fig. da) an d then for that shown 
in Fig. IHb). A positive infinitesimal V in Eq. (jUJ) is replaced by T = 0.01J in the ci-wave 
pairing state and T = 0.05 J in the normal state. Although the choice of a finite T is done for 
numerical convenience, it may simulate damping of quasiparticles by static defects in real 
materials and broadening due to limited energy resolution in inelastic neutron scattering 
experiments. 

Figure El shows (q, uj) maps of the imaginary part of x(q, tu) for the FS shown in 
Fig. E^a) for several choices of temperatures; q is scanned along (0.47T, vr) < q < (tt, tt) and 
(it, tt) > q > (0.57T, 0.57r); along each direction the highest weight positions of Inr\;(q, oj) are 
represented by cross symbols; in Fig. Of a) we also plot the threshold energy of individual 
particle- hole excitations by a gray line; Figures HHa) and (b) are composed of three different 
energy regions with an optimal color map scale. There appear gapless excitations along 
q = (q,q)/y/2 in Fig. |3fa). This diagonal IC signal is due to scattering processes between 



28, 



With 



gap nodes of the ci-wave singlet pairing, as already clarified about a decade ago. 
increasing u, scattering processes around IC positions at q = (tt ± 2irr], tt) and (tt, tt ± 2ttt]) 
begin to contribute. The highest spectral weight positions (cross symbols) appear close to 
the threshold energy for both the IC and the diagonal IC signals; the IC signal becomes 
stronger than the diagonal IC signal. The incommensurability 77 typically tends to decreases 
with uj, resulting in a commensurate peak at a particular energy (a; ~ 0.45 J), which we de- 
fine as cuq] note that ujq corresponds also to a peak energy of Imx(Q, uj) in the <i-wave pairing 
state. Looking closely at a region near ujq, we find that the peak position of Imx(q, uj) is 
located slightly below the threshold energy, where the denominator in Eq. (JHJ) vanishes, that 
is, Imx(q, uj) has a pole. Hence the strong signals near uj ~ u>q are interpreted as collective 
particle-hole excitations, namely the so-called resonance. But the peak position appears 
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very close to the threshold energy, not well inside the gap. Thermal broadening and the 
spectral weight broadening due to a finite T [Eq. (JHJ)] easily make the resonance overdamped 
through mixing with individual particle-hole excitations. In addition, the (q, uj) region where 
the resonance is realized is very limited. Therefore, particle-hole excitations below ujq are 
mainly characterized by a downward dispersion of individual excitations, which traces the 
q dependent threshold energy. 

Above ujq, the strong spectral weight appears around q = Q (cross symbols). We can, 
however, read off from Fig. a) that such a commensurate signal is accompanied by IC and 
diagonal IC substructures for uj > 0.5 J, which disperse outwards with increasing uj. 

The (q, uj) map in Fig. HD^a) is robust against temperature as seen in Fig. [31(b) and 
the downward dispersion below ujq remains even in the normal state, namely above 
7rvb(= 0.104J) [Fig. Cfl[c)]; Trvb is the onset temperature of the o?-wave pairing gap, and 
is interpreted as the pseudogap crossover temperature T* in the underdoped regime and 
as the superconducting phase transition temperature T c in the overdoped regime of high-T c 
cuprates. [5_4j Although IC and diagonal IC peaks are seen in Fig. Efc), the peak structures 
are very broad because of thermal broadening. When temperature is reduced under the 
condition A = (x T and /i are determined self-consistently) , we see clearer IC and diago- 
nal IC signals in Fig. Efd); the former is stronger than the latter. These incommensurate 
structures in the normal state can be traced back to the FS geometry shown in Fig. ^a). 
Besides the presence of incommensurate nesting vectors, there are no particle-hole scattering 
processes across the FS with q = Q and a small uj, which substantially reduces the spectral 
weight around Q, leading to robust incommensurate structures of magnetic excitations. In 
Figs- Etc) and (d), we see a IC substructure above uj ~ 0.2J along q = (q x ,7i), which dis- 
perses outwards with increasing uj. This substructure becomes stronger for a larger 5 and a 
smaller V. 

Next we analyze the S dependence of Imx(q, uj). Figure |U shows (q, uj) maps of Inr\;(q, uj) 
in the ci-wave pairing state (T = 0.01J). A downward dispersive feature below ujq still 
appears for a lower 5, accompanied by a reduction of ujq. Above ujq, on the other hand, 
an upward dispersive feature shows up with decreasing 5. In particular, its spectral weight 
becomes larger than that of the downward dispersion in Fig. Efb) (see values of color map 
index). This upward dispersion is not due to collective excitations, but just a peak of 
individual particle-hole excitations. In the vicinity of the AF instability (5 ~ 0.02), however, 



S 



this upward dispersion changes into an overdamped collective mode in the sense that the real 
part of the denominator of Eq. (jHJ) vanishes. In the normal state, low-energy incommensurate 
structures become less clear for a lower 5 since the FS geometry tends to allow particle-hole 
scattering processes with q = Q on the FS. 

Results similar to Figs. 01 and 0] are obtained for the FS shown in Fig. ^b). (q, uj) maps 
of Imx(q, uj) are shown in Fig. El for several choices of T. In Fig. Efa), peak positions of 
Imx(q, uj) (cross symbols) trace a downward dispersion below ujq(= 0.49J), which comes 
from individual particle-hole excitations. While a collective mode, namely the resonance, 
does not appear in Fig. Efa) (5 = 0.15), it can appear for a lower doping region. But even 
in this case, we checked that the collective mode is not well separated from the continuum 
excitations and thus easily overdamped for a finite T and T. Moreover the collective mode 
appears in a very limited (q, uj) region as seen in Fig. [31(a) . In this sense, the downward 
dispersion below ujq is characterized mainly by individual excitations even for a lower doping 
rate. Although spectral weight above ujq is characterized by a commensurate signal at 
q = Q, we see a small segment of an upward dispersion just above ujq in Fig. Efa). When 
the hole density is decreased, this upward dispersion becomes clearer and a result similar to 
Fig. 0] is obtained. Overall features seen in Fig. OJa) still survive for a higher temperature 
[Fig. 03(b)], although the temperature in Fig. |3fb) is high enough to smear the segment of an 
upward dispersion for uj > ujq. Even above T RV b(= 0.104 J), incommensurate structures are 
still seen in Fig.E^c) and become clearer when temperature is reduced under the condition of 
A = [Fig.Efd)]. The robust property of the incommensurate structures in the normal state 
[Figs. Efc) and (d)] can be understood in terms of the geometry of each deformed FS shown 
in Fig. ^b) by the same argument as that already given in the context of Fig. El In reality, 
there is a weak interlayer coupling and thus the resulting FSs are composed of a holelike FS 
and an electronlike FS as shown in Fig. |21 which then allows low-energy scattering processes 
with q = Q through the interlayer coupling. Such effects are, however, too weak to smear 
incommensurate structures for a realistic parameter as was explicitely calculated in Ref . j3| . 
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IV. DISCUSSION 



A. Magnetic excitations in La-based cuprates 

Now we compare our results with neutron scattering data. For both FSs shown in 
Figs, IHa) and (b), we have obtained IC magnetic excitations for a relatively low uj not 
only in the <i-wave pairing state but also in the normal state [Figs. H3 and Ej, which well cap- 



DDDD 



tures the most prominent feature of magnetic excitations in La-based cuprates. j 

The incommensurability typically decreases with increasing u and the signals finally 
merge into a commensurate peak at u — uq. This dispersive feature is consistent with 
recent data for LBCO with x = 0.125,0 and also for LSCO with x = 0.10 and 0.160 
where measurements were performed up to a certain energy below ujq. 

For a higher energy region in LBCO with x = 0.125, Tranquada et a/.^| reported 
upward dispersive signals of Imx(q, oS). Their obtained spectral weight distribution was 
very broad with substantial spectral weight at q = Q, and the intensity difference between 
IC (or diagonal IC) positions and a commensurate position was very small. Our results for 
5 = 0.15 (Figs. |3]and|5J) show typically a broad spectral weight around q = Q for a high 
energy region, but an upward dispersive feature shows up with decreasing 5 in the ci-wave 
pairing state (Fig.©. Although their measurementQ was performed above T c , we may 
assume the data were obtained in the pseudogap state in LBCO, which is associated with 
the d-wave pairing state in the slave-boson mean-field theory. 0] Under this assumption, we 
can capture their data within the present study. 

Because of rfFSD correlations, the shape of the FS in La-based cuprates may depend 
strongly on the crystal structure. 0] In the presence of the LTT structure, which was ob- 
served in LBCO0, Nd-doped LSCO 6^, and LSCO with x = 0.120, we expect that 
rfFSD correlations lead to a strongly deformed FS as shown in Fig. El through coupling to 
a small xy anisotropy of the lattice. The corresponding magnetic excitations are shown in 
Fig. El Even in the absence of the (static) LTT, the soft phonon mode toward the LTT 
structural phase transition, whose energy u p h is about a few meV, was observed in LSCO 
with x < 0.180 HH0 and LBCO with x = 0.125 79]. In this case we expect dy- 
namical fluctuations of the <iFSD within a time scale shorter than uj~^. Thus high energy 
probes such as ARPES and neutron scattering may see an instantaneous rfFSD. On the 
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other hand, in the absence of the LTT and its slow fluctuations, which we expect roughly 



for T > 100 - 200K, 



H 76, 77,0 



791 ] no driving force leading to a strongly deformed FS 



may be present. We expect a FS as in Fig. [T] (a) and magnetic excitations as in Figs. El 
and I3J In the present ciFSD correlation scenario, therefore, the change of the FS shape is 
predicted as function of T. The ARPES data for high T have not yet been obtained for 
La-based cuprates and this prediction can be tested in a future. 

An indirect evidence of the FS change was recently obtained in a neutron scattering 
experiment for LBCO with x = 0.125. This material shows almost a first-order-like LTT 
structural phase transition at 60K.jzj| Fujita et al^ measured a temperature dependence 
of the incommensurability 77 at a low energy and found a sizable change of rj at the LTT 
transition. This is naturally understood in terms of the (static) FS change scenario from 



Fig. to Fig. W[b) below 60K; the FS shown in Fig. Cfb) tends to favor a larger rj, the 
same tendency as the experiment. The authors in Ref. 9|, however, interpreted the data 
differently as a lock-in effect that the periodicity of charge stripes tends to be commensurate 
with the lattice potential of the LTT structure, although a charge order signal developed 
gradually below 50K, not directly below the LTT transition temperature. 

The state with a ci-wave deformed FS has the same symmetry as the so-called electronic 
nematic state. 8(| The nematic order was extensively discussed for cuprates. j3| However 
this nematic order was discussed as coming from partial spin-charge stripe order, not from 
(iFSD correlations; magnetic excitations were then discussed in terms of stripes. Our rfFSD 
does not require charge stripes nor their fluctuations, but is driven by forward scattering 
processes of quasiparticles, which provides another route to the nematic state. An interesting 
open question is whether the e?FSD state might have an instability toward a charge ordered 
state such as stripes. Even if it were the case, our calculation shows that many prominent 

1-captured 



2J, I25J do 



features of magnetic excitations observed in La-based cuprates are already we 
without stripes. Observations of weak charge order signalsjj], 0, 21, 0, 0, 
not necessarily mean that charge stripes are crucial to magnetic excitations. Effects of a 
charge order on magnetic excitations can be higher order corrections beyond the present 



renormalized RPA and might be responsible for a realization o 
order around hole density 1/8 in La-based cuprates, 
captured in the present framework. 



21 



static IC antiferromagnetic 



24, 



81 



which cannot be 
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B. Comparison with magnetic excitations in YBCO 



We finally discuss how similarities and differences of magnetic excitations between La- 
based cuprates and YBCO arise from the present framework. In the slave-boson mean-field 
theory, material dependences are described by different choices of band parameters. For 
La-based cuprates, we have taken t/J = 4, t'/t = —1/6, t"/t = 0, and r = 0.35 (see Sec. II) 
while t/J = 2.5, t'/t = -0.30, t"/t = 0.15, and r = 0.5 were invoked for YBCO (see 
Ref . [44] , where the bilayer coupling was also taken into account and a comprehensive study 
of magnetic excitations in YBCO was performed including effects of the (iFSD). 

A well-known distinction between La-based cuprates and YBCO is that low-energy IC 



magnetic signals are realized up to a temperature much higher than T c in the former 

| while they are realized only below T c or possibly below the pseudogap temperature T* 



in the latter 



aey are i 

•BO, El 



This difference comes from a FS difference due to different choices of 



t'/t and t"/t. For the FS shown in Fig. Q there are no particle- hole scattering processes 
with q = Q for low energy, yielding robust incommensurate structures even in the normal 
state, while the FS for YBCoQ allows such scattering processes, smearing incommensurate 



28. 20|, to understand this 



features in the normal state. We confirm early works[27, 
material dependence of magnetic excitations. 

In the d-wave superconducting state, on the other hand, both LSCO 



and YBCO 



5, 



14l Il5j show IC magnetic signals for low uj and its incommensurability rj decreases with 



increasing to. This behavior is well-reproduced in the present framework, and does not 
depend on the choice of band parameters as emphasized before. 

The so-called resonance sig nals of magnetic excitations were reported for YBCO in the 



superconducting state 



10 



11 



12] and a pseudogap stateQ, I>u1 uol for LSCO. In the 



present theory, this resonance is interpreted as collective particle-hole excitation as a. 



discussed by many authors. [ 



40 



41 



42, 



43, 



ready 



4J The 



0, 32, 0, y, y, 3 y y , y, 

realization of collective excitations depends strongly on the choice of r in Eq. (JHJ); it is more 
favorable for a larger r. Although our relatively small r for LSCO can produce a resonance, 
its energy appears very close to the threshold energy of continuum excitations, making the 
resonance easily overdamped by thermal broadening and spectral broadening due to a finite 
T [Eq. (JHJ)]. Therefore the resonance signals are not so clear for LSCO as for YBCO. It is 
to be noted that if we invoke a larger r for LSCO, keeping other parameters unchanged, 



12 



well-defined resonance signals appear as in the case of the parameters for YBCO.j^J There 
is another factor, a FS difference, which is less effective than a factor r. Since the FS 
invoked for LSCO (Fig. is closer to the (tt/2,tt/2) point than that for YBCO, lowest- 
energy particle-hole scattering processes with q = Q appear at k x = k y when a large <i-wave 
paring gap opens on the FS. For such scattering processes, the coherence factor C* kk+q 
[Eq. (0)] vanishes. Hence at zero temperature the imaginary part of x (Q?^) does not 
show a jump at the threshold energy, but increases continuously; the corresponding real 
part of x (Q,c<j) does not have a log- divergence there. This is the case for 5 < 0.08 for 
the prese nt p arameters for La-based cuprates. Since one of well-known mechanisms of the 



resonance 



it paramcti 
H 32, 3, 



34. 



35 



36 



37 



39 



40 



41 



42 



43 



441 ] is based on this log- divergence 



of Rexo(Q)^)) the resonance is further disfavored for La-based cuprates. 

Above the resonance energy, the so-called upward dispersive mode was found in YBCO 



15l . \\\ , which was captured in the present slave-boson mean- field framework 4^] as well as 







other phenomenological calculations. 



35, 



A similar high energy spectrum, which was 
much broader than that in YBCO, was also reported for LBCO with x = 0.125. ^| Our 
results in Figs. 13141 and El capture this experimental data. The realization of an upward 
dispersion becomes more favorable for a smaller 5, a larger r, and a smaller t/J. Since 
the latter two factors are not respected for the present parameters for La-based cuprates, 
the origin of the upward dispersion is different from that for YBCO: a peak of individual 
particle-hole excitations (except for the vicinity of the AF instability) for La-based cuprates 
and overdamped collective excitations j^J for YBCO. 



V. CONCLUSION 



We have computed Imx(q, uo) up to a high energy region (< J) and shown that salient 
features of magnetic excitations in La-based cuprates are well-captured in terms of particle- 
hole excitations. IC magnetic signals are realized for low u in both the (i-wave pairing state 
and the normal state, and they disperse with increasing u, merging into a commensurate 
signal at u — cjq. The resonance mode can be realized close to u — o>q, which is easily 
overdamped by mixing with individual particle- hole excitations for finite T and T. Above o>q, 
upward dispersive features of Inrx(q, u) appear especially for a low 5 in the <i-wave pairing 
state, while typically broad spectral weight distribution is obtained around q = Q in the 



13 



3,0 



normal state. In the present theory, <iFSD correlations |48l |60( have played an important 
role for the shape of the FS. In particular, we expect a change of a (static) FS shape 
across the LTT structural phase transition, which in general yields a sizable change of 
incommensurability. Such a change was recently observed in LBCO. 9] Combined with 
previous work 3| focusing on the low energy region of magnetic excitations, the present 
study confirms the crucial role of particle-hole excitations up to high energy for magnetic 
excitations in La-based cuprates. 
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FIG. 1: Fermi surface for 5 = 0.15: (a) isotropic case and (b) anisotropic case with a = 0.85, where 
two Fermi surfaces (solid line and gray line) are superimposed. The Fermi surface is defined by 
Ck = at T = 0.01J in the d- wave pairing state. 




FIG. 2: Alternate stacking of ci-wave deformed Fermi surfaces along the z axis in the presence 
of the LTT in La-based cuprates. Because of a weak interlayer coupling t z , the resulting Fermi 
surfaces (right panel) are composed of a holelike Fermi surface (solid line) and an electronlike Fermi 
surface (gray line). 



19 




(0.4jt,jt) (Q.5jt,0.5rt) (0.4jt,JT) (0.5jt,0.5jt) 




t (?»,3I) 0?,g)/V2 f f ( 9i .,3t) (31, It) (*g)/V2 f 

(0.4JT,JT) (0.5jt,0.5jt) (0.4lt,Jt) (0.5JI, 0.5jt) 



FIG. 3: (color online) (q, uj) maps of Imx(q, u) for 5 = 0.15 at several choices of T for the Fermi 
surface shown in Fig. ^a); Tryb = 0.104 J, but A = is assumed in (d); the q scan directions are 
(0.47T, vr) < q < (jr, tt) and (tt, tt) > q > (0.57T, 0.57r); the cross symbols represent the highest weight 
positions along q = (q x ,Tv) and (q,q)/y/2; the gray line in (a) is a lower edge of the continuum 
excitations; (a) and (b) are made of three different energy regions with an optimal color map scale. 
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FIG. 4: (color online) (q, u) maps of Imx(q, uS) at T = 0.01J for 5 = 0.10 (a) and 0.05(b) 
corresponding result for 5 = 0.15 is shown in Fig. |2fa). 
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FIG. 5: (color online) (q, uj) maps of Imx(q,Lj) for 5 = 0.15 at several choices of T for the Fermi 
surface shown in Fig. ^b); Trvb = 0.104 J, but A = is assumed in (d); the corresponding results 
for the Fermi surface in Fig. ^a) are shown in Fig. |21 
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